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Abstract 

We perform a brief but critical review of the Landauer picture of transport that clar- 
ifies how decoherence appears in this approach. On this basis, we present different 
models that allow the study of the coherent and decoherent effects of the interaction 
with the environment in the electronic transport. These models are particularly well 
suited for the analysis of transport in molecular wires. The effects of decoherence 
are described through the D'Amato-Pastawski model that is explained in detail. We 
also consider the formation of polarons in some models for the electron- vibrational 
interaction. Our quantum coherent framework allows us to study many-body inter- 
ference effects. Particular emphasis is given to the occurrence of anti-resonances as 
a result of these interferences. By studying the phase fluctuations in these soluble 
models we are able to identify inelastic and decoherence effects. A brief description 
of a general formulation for the consideration of time-dependent transport is also 
presented. 
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1 Introduction. 



Electronic transport in biological and organic molecules has become a very 
exciting field [1,2] that brings together ideas and results developed during the 
past two decades in many branches of Physics, Chemistry and Biology. On the 
basis of this synergic interaction one can foresee very innovative results. Much 
of its wealth comes from the fact that molecules are intrinsically quantum 
objects, and quantum mechanics always defies our classical intuition. When 
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this happens, we are almost certain to find new "unexpected phenomena" 
leading to prospective applications. In turn, with a few exceptions, it is not 
yet clear how Nature exploits quantum efi:ects in biological systems. However, 
it seems hkely that evolution has made use of the details of the quantum 
tunneling process in modulating charge transport [3]. Hence, by pursuing an 
exploration of dynamical quantum phenomena at the molecular level one may 
also expect to be more prepared to discover Her hidden ways. 

A paradigmatic example of the above situation is the electron transfer in DNA 
strands [2] . There is evidence that at least two mechanisms are present in this 
case: 1) A coherent tunneling between the base pairs constituting the donor 
and acceptor centers. 2) An inelastic sequential hopping through bridging cen- 
ters. Theoretically, there is a need to unify the description of these extreme 
regimes as well as to study the possible role of vibrations and distortions of the 
DNA structure. Our present understanding of these basic physical processes 
comes from the field of Quantum Transport, which evolved from the descrip- 
tion of non-crystalline materials [4] and reached its climax with nanoelectronic 
[5] devices. Since these last structures may have a very small length scale the 
term "artificial atoms" [6] is amply justified. Many of the phenomena appear- 
ing in these systems could be summarized in the fact that the wave nature 
of electrons made them susceptible to interference phenomena. These can be 
assimilated to the propagation of light in some complex Fabry-Perrot inter- 
ferometer producing either well defined fringes or "speckle-patterns" . Failures 
in that simplistic description is essentially due to interactions with excessively 
complex environmental degrees of freedom (such as thermal vibrations), and 
our lack of control over them is interpreted as "decoherence" [7] . This often 
justifies the use of a classical description of the transport process. However, 
if one is to borrow some idea from the theory of solid state physics, it could 
be that coherent interactions between electrons and lattice distortions give 
rise to new effects such as assisted tunneling [8] or even superconductivity. 
Indeed both phenomena could give striking results in organic systems [9,10]. 
One is then compelled to develop new tools to describe electron-phonon (e-ph) 
coherent processes in molecular devices. 

In what follows we will make a brief description of quantum interference phe- 
nomena establishing a ground level language based on the simplest physical 
models. While following our personal pathway of many years through the gen- 
eral ideas of transport in the quantum regime, we expect to induce a new 
perspective into the subtle mechanisms that lead to the degradation of the 
simple quantum effects through the interactions with the environment (i.e. 
decoherence). As a token, we will visualize coherent effects emerging from the 
electron-lattice interaction that can be exploited in new useful ways. 

In Section 1 we recall the Landauer ideas [11] for transport which will also 
serve to adopt a basic language and give a conceptual framework into which 
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the phenomenological aspects of decoherence can fit. In Section 2 we review 
the D'Amato-Pastawski (DP) model [12] for coherent and decoherent trans- 
port commenting its strengths and hmitations. In Section 3 we introduce a 
polaronic model that hints at the properties of the complete electron-phonon 
Fock space. This not only sheds light on the decoherence process but can also 
be used to predict new phenomena such as the coherent emission of phonons. 
We devote Section 4 to a brief description of the time dependent transport. 
Section 5 gives a general perspective. 



2 Landauer's picture for transport. 



The framework that inspired most of the described experimental developments 
in mesoscopic electronics was the Landauer's [11] simple but conceptually new 
approach. Besides the "sample" or device, he explicitly incorporated the de- 
scription of the electric reservoirs connected to "the sample" . The role of reser- 
voirs can be played not only by electrodes but also could be spatial regions 
(localized LCAO) where the electrons lose their quantum coherence. This last 
situation describes transport in the hopping regime [14] (we will see more on 
this latter). The simplest mathematical description is obtained if one thinks 
them as one- dimensional wires that connect the individual orbitals i to elec- 
tronic reservoirs characterized by the statistical distribution function U{e). In 
that case, the electronic states are plane waves describing the different bound- 
ary conditions of electrons "in" or "out" the reservoir. An electron "out" from 
reservoir connected to "site" i has a probability Tj^i{e) to enter the reservoir 
connected to "site" j. A representation of such a situation for the case of 
three reservoirs is sketched in Fig. 1. The current per spin state at reservoir j 
is obtained by the application of the Kirschhoff 's law (i.e. a balance equation 
[13]): 



The meaning of this equation is obvious: it balances currents. Each reservoir j 
emits electrons with an energy availability controlled by a Fermi distribution 
function ii^e) = l/[exp[(£ — fii)/kBT + 1], with a chemical potential, Hi ~ 
/lo+Sui, displaced from its equilibrium value Ho- One can assimilate Vj = Sni/e 
as voltages. The density of those outgoing states is |iVj(e) (half the total) and 
their velocity Vi. It was essential in Landauer's reasoning to note that in a 
propagating channel the density of states Ni is inversely proportional to the 
corresponding group velocity: 

Ni = l/{vih). (2) 
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Fig. 1. Representation of a three probe measurement. The voltmeter may be strongly 
coupled and is a source of decoherence. 

This is immediately satisfied in one-dimensional wires but its validity is much 
more general. This fundamental fact remained unnoticed in the early discus- 
sions of quantum tunneling [5] and it is the key to understand conductance 
quantization. The Gj^j = ^(g) are hence Landauer's conductances per spin 
channel. For perfect transmitting samples ^ is either 1 or and one obtains 
the conductance quantization in integer multiples of /h. 

Notice that there is no need for the traditional [1 — ij{e)] factor to exclude 
transitions to already occupied final states. In a scattering formulation, any 
"in" state contains a linear combination of "out" states. Although two differ- 
ent "in" states (e.g. on the left and right electrodes) could end in the same 
final "out" state, unitarity of quantum mechanics assures that both sets are 
orthogonal. Here the transmission coefficients may depend on the external 
parameters such as voltages. 

An important particularity of our Eq. (1) is that it does not exclude sites i = j 
from the sum. This contrasts with the original multichannel description [13] 
and is of utmost importance in the treatment of time dependent problems 
where Vj = Vj(t) as will be seen in Section 4. Our Eq. (1), always used with 
Eq. (2), has a full quantum foundation within the Keldysh formulation of 
Quantum Mechanics (e.g. see Eqs. (5.6-7) in ref. [15]). It can be fully expressed 
in terms of quantities obtained from a Hamiltonian model such as local density 
of states and Green's functions. We adopt here a notation consistent with these 
formal developments. 

2.1 Phenomenology of decoherence. 

A first alternative to include decoherence in steady state quantum transport 
was inspired in the Landauer's formulation. There, the leads, while accepting 
a quantum description of their spectrum of propagating excitations, are the 
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ultimate source of irreversibility and decoherence: electrons leaving the elec- 
trodes toward the sample are completely incoherent with the electrons coming 
from the other electrodes. In fact, it is obvious that a wire connected to a volt- 
meter, by "measuring" the number of electrons in it, must produce some form 
of "collapse" of the wave function leading to decoherence (see Fig. 1). Be- 
sides, no net current flows toward a voltmeter. The leads are then a natural 
source of decoherence which can be readily described in the Landauer picture 
if one uses the Landauer conductances together with the Kirschhoff balance 
equations. This fact was flrstly realized by M. Biittiker [16]. Let us see how 
it works for the case of a single voltmeter in the linear response regime. In 
matrix form: 



II 



\}r) 



— [Gr,l + G<^,l] ^\.,4> Gl,r 

G</,,L — [Gr,</, -I- Gl,,/,] G<^,r 
Gr,L Gr^0 — [G03 -I- Gl,r] 



(3) 



Here the unknowns arc II, Ir and V<^ = The second equation must be 

solved with the voltmeter condition L = 0: 



6 6 

= j^T^,L(SiJi^ - 5 III,) + -Ti^^^{5iJL^ - Siir), 



(4) 



giving us 6 11(1, to be introduced in the third equation: 



6 6 

Ir = ^7r,l(5//l - Shr) + -Tn,4,{6iJ,4, - 5/xr) (5) 



to obtain the current 



with 



6 ~ 



^R,(a7'(A,L 



-Lr,l — -L R,L + 7f —pf — 

The first term can be identified with the coherent part, while the second is the 

incoherent or sequential part, i.e. the contribution to the current originated 
from particles that interact with the voltmeter. This corresponds to an effective 
conductance of 

Gr,l = Gr,l + (Gi,,^-^ + G^,l-')-\ (7) 

which can be identified with the electrical circuit of Fig. 2. This classical view 
clarifies the competition between coherent and incoherent transport. However, 
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Fig. 2. Classical circuit representation of the non-classical system of Fig. 1. 

this circuit docs not imply that in Quantum Mechanics one cannot modify one 
of the resistances without deeply altering the others. This fact can become very 
relevant in partially coherent regimes. 

So far with the phenomenology. The next important step is to connect these 
quantities with actual model Hamiltonians. This connection was made explicit 
by the contribution of D'Amato and Pastawski [12]. 



3 The D'Amato-Pastawski model for decoherence. 

3.1 Reducing the Hamiltonian 

Let us first review the basic mathematical background that made possible the 
selection of a simple Hamiltonian that best represents the complex sample- 
environment system. The objective was to find a simple way to account for the 
infinite degrees of freedom of a thermal bath and/or electrodes and use the 
exact solution in the Landauer's transport equation. First, we recall that one 
can always eliminate the microscopic degrees of freedom [17,18] generating an 
effective Hamiltonian, expressed in a basis of localized orbitals. This produces 
effective interactions and energy renormalizations which depend themselves on 
the observed energy. Furthermore, one can include a whole lead in a Hamil- 
tonian description through a correction to the eigen-energies which has an 
imaginary part. In fact, an electron originally localized in the region called 
"the sample" should eventually escape or decay toward the electrodes. In a 
microscopic description this is equivalent to the decay of an excited atom ac- 
cording to the Fermi Golden Rule. Hence, there is a escape velocity associated 
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with the energy uncertainty of a local state: 

where a is a lattice constant. Let us see a simple example. Consider a sin- 
gle LCAO which could be the highest occupied molecular orbital (HOMO) 
or the lowest unoccupied molecular orbital (LUMO) responsible for resonant 
transport 

coupled with a qyia.si- continuum of electronic states of the electrode. They are 
described by a "lead" connected to the voltmeter: 

k 

The sample-lead interaction term is: 

"^O-lead = J2 {^kfiCkCo + C.C.) . 
k 

The unperturbed atomic energy Eq will become corrected by the presence of 
the lead. In second order of perturbation we get 

i-'^Eo^(^) = hm^ Y: = '""^^^oiEo) T i"=^'^ro(i?o). (9) 

The sign -|- or — of the infinitesimal imaginary energy 77 is introduced to 
handle eventual divergencies in the sum. Since it determines the sign of time 
in the evolution, the supra-indcx or ^4 corresponds to either a retarded or 
advanced propagation. The imaginary component appears because the elec- 
trode spectrum, described by the density of states A^icad(£^)j is continuum in the 
neighborhood of energy Eo + /S.{Eo). This makes possible the irreversible decay 
into the continuum set of states k, not included in the bounded description. 
The general functional dependence of Ao(£^o) ~ iro(-E'o) is better expressed in 
terms of the Wigner-Brillouin perturbation theory: 

/oo I T/, I ^ 
l^N,,^{Ek)dEk, (10) 
-00 £ — -C/fe 

where p stands for principal value and iViead(-E'fc) is the density of states at the 
lead. Similarly: 

/°° 
\Vkfl\ Nie^A{Ek)5[e-Ek]dEk. (11) 
-00 

The evaluation of Eq. (11) and (8) at £ = £^0 constitutes the Fermi Golden 
Rule (FGR). Of course these quantities satisfy the Kramers- Kroning relations 
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Me) - Ipf "^ie'. (12) 

An explicit functional dependence on the variable £ contains eventual non- 
FGR behavior. The FGR describes reasonably electrons in an atom decaying 
into the continuum (electrode states), or propagating electrons decaying into 
different momentum states by collision with impurities or interaction with a 

field of phonons or photons. In some of these cases, we have to add some 
degrees of freedom to the sum (the phonon or photon coordinates). A process 
a may produce contributions "S^ to the total self-energy S^. 

^Ointeraction^'^O = + (13) 

with 



Eo^(^)(£) = (Ao(£)Tiro(£))4?o 

= ^("Ao(£)Ti"ro(£))c+Co. 



(14) 
(15) 



Besides, the best way to do perturbation theory to infinite order is the frame- 
work of Green's functions. 

The unperturbed retarded Green's functions are defined as the matrix ele- 
ments of the resolvent operator 



G°^{£) = lim [(£ + ir])i - HT^ 
»?— >o+ 



(16) 



for the isolated sample. It is practical to represent this as a matrix, G°^{e), 
whose elements are written in terms of the eigen-energies and eigen- 
functions tjj^^^ (r) = I]j Ma,i97i(r) of the isolated sample as 



G^f (£) = lim ^ ^ 



a 



'?^o+ a £ + - ^; 

The Local Density of States at orbital i is calculated as 



(17) 



(18) 



The Fourier transform 



(19) 
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Fig. 3. Feynman representation of the Dyson equation for the single particle Green's 
function (line). An electron-phonon self-energy is evaluated in terms of the sample 
electron (line) and phonon (wave) Green's functions. The leads self-energies contain 
the hopping (dot) and the propagator in the lead. 



is solution of: 

d 



SiAt2-ti) (20) 



with ^2 > ^1- The identity matrix is = Sij. Hence G°j{t2 — ti) is — times 
the probability amplitude that a particle placed in j-th orbital at time ti be 
found in i-th one at time t2- The advantage of the method is that Gf^it) can 
be calculated numerically from Gfj{e) without a detailed knowledge of the 
eigen-solutions of the perturbed system. For G^j{e) one uses: 

G^(£) = [£l-(H„ + E^(£))]-^ (21) 

oo 

= ^^(e) + G^'^ie) J2 [S^ (6)0"^ (£)]"" (22) 

n=l 

= G<'-^(£) + ^^(s) (£)G^(£) (23) 

In the second line we have written the usual forms of the Dyson equation from 
which the Wigner-Brillouin perturbative series and its representation in Feyn- 
man diagrams can be obtained. In Fig. 3 we show the graphical representation 
of the Dyson equation considering two local contributions to the self-energy: 
the electron-phonon interactions and the escape to the leads. 

For a brief tutorial on the calculation of the Green's function in discrete sys- 
tems see Ref. [19]. Although the formalism seems to introduce some extra 
notation it has various conceptual advantages. For example, it is straightfor- 
ward to use Eq. (21) to prove the optical theorem [19]: 

[G^ - G^] = G^ [E« - E^] G^ (24) 

of deep physical significance since it is an integral equation relating the local 
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densities of states given by Eq. (18) and the decay rates provided by Eq. (14). 
Perhaps the most important advantage, is that they can be used also in the 
Quantum Field Theory [20,21] to deal with the many-body case. 



The direct connection between observables and the transmittances used in 
the Landauer formulation and the Green's function follows intuitively from the 
Green's function probabilistic interpretation. This was formahzed by D'Amato 
and Pastawski who used a result of Fisher and Lee [22] that related the trans- 
mittance to a product of the Green's function connecting two regions and 
their group velocities. By noting the discussed relation between the group ve- 
locity and the imaginary part of the effective potential established in Eq. (8), 
DP obtained a simple expression for the transmittance, which in our present 
notation can be written as 

Tja,(^{e) ^ = [2"r,(£)] Gf/e) [2%{e)] Gf^^ie)] . (25) 

The left supra-index a in °T indicates each of the independent processes pro- 
ducing the decay from the LCAO modes (right subindex) associated with the 
physical channel. In the basis of independent channels the complex part of the 
self-energy is diagonal. Hence the sum over initial states and processes in the 
physical channels R and L of Eq. (1) determines a conductance: 

Gr,l = ^4TV[rR(6) Gl^ie) rUe)Gl^{e)]. (26) 

This matrix expression is simply a compact way to write Eq.(l), using Eqs. 
(25) and (8) together to compute the linear response conductance between any 
pair L and R of electrodes. The sum of initial states at left (L) is the result 
of the product of the diagonal rL(^) form of the broadening matrix, while the 
final trace is the sum over final states at right (R). With different notations, 
these basic ideas recognize by now many applications to the field of molecular 
electron transfer [23] . 

Originally, Fisher and Lee considered only the escape velocity to the leads (i.e. 
a = f3 = lead). Our point is that any other process which contributes to the 
decay giving an imaginary contribution to the self- energy would be described 
by Eq. (25). In particular, this will be true for a "decoherence" velocity that 
degrades the coherent current. This view was proposed in DP [12] adopting a 
discrete (tight-binding) description of the spatial variables at each point (or 
orbital) in the real space a decay rate was assigned which is balanced by the 
particles reinjection in the same site, described as local reservoirs. 

We now show how the DP model applies to a simple tunneling system. Con- 
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sider the unperturbed Hamiltonian, 



Tio = E + E [n^o + . (27) 

Since the indices i and j (sites) refer to any set of atomic orbitals, the inter- 
actions are not restricted to nearest neighbors. However, for the usual short 
range interactions, the Hamiltonian matrix has the advantage of being sparse. 
The local dephasing field is represented by 

N 

't't^ = ^ _i c+q (28) 
1=1 

where '^V — h/{2T^). We consider for simplicity only two one-dimensional 
current leads L and R connected to the 1st. and the A'"th orbital states respec- 
tively, 

leadsj^R ^ _• j^Lp p g+g^j _ ^29) 

We see that the 1-st state has escape contributions both, toward the current 
lead at the left, ^Fi, and to the inelastic channel associated to this site, "^Fi. 
The on-site chemical potential will ensure that no net current flows through 
this channel. 



3.2 The solution for incoherent transport 



To simplify the notation we define the total transmission from each site as: 



N+1 

E Tj, 

j=0 



AttNi^Ti fori = 

47rAr>Fi for 1 < i < AT 
47rA^jv^rAr fori = A^+l 



(30) 



The last equality follows from the optical theorem of Eq.(24). The balance 
equation becomes 



N+1 



Ij = = - (l/gi) S/jLi + E TijS/jLj, 



(31) 



j=0 



where the sum adds all the electrons that emerge from a last collision at other 
sites (j 's) and propagate coherently to site i where they suffer a dephas- 
ing collision. These include the electrons coming from the current source i.e. 
Ti^L^fJ^h and the current drain. However, since we refer all voltages to the last 



11 



one, Tj3(5//R = 0. The first term accounts for all the electrons that emerge 
from this collision on site i to have a further dephasing collision either in the 
sample or in the leads. The net current is identically zero at any dephasing 
channel (lead) . The other two equations are 

TV 

II = -I = - (1/^l) + E TlJ^j, (32) 

j=0 
N+l 

Ir = I = - (I/s'r) ^A^r + E TR,j6fXj. 

j=0 

Here we need the local chemical potentials which can be obtained from Eq. 
(31). In a compact notation, these coefficients can be arranged in a matrix 
form which excludes the leads that are current source and sink: 

1/ 9i ~ ^1,1 ^1,2 2 ■ ■ ■ Ti TV 

W = T3,i n,2 1/93 - 7^3,3 • • • n,N ■ (33) 

Tn,1 Tj^,2 Tj^,3 • ■ ■ ^/ 9n — Tn,n 

Then, the chemical potential in each site can be calculated as 



N 



3=^ 



(34) 



Replacing these chemical potentials back in Eq. (31) the effective transmission 
can be calculated 



TV AT 

7R,L = Tk,L + EE^R,.[w-' 

j=l i=l 



T, 



3,t 



i,L- 



(35) 



The first contribution in the RHS comes from electrons that propagate quan- 
tum coherently through the sample. The second term contains the incoherent 
contributions due to electrons that suffer their first dephasing collision at site 
i and their last one at site j. 

Until now the procedure has been completely general, there is no assumption 
involving the dimensionality or geometry of the sample. The system of Fig. 4 
was adopted in DP only because it has a simple analytical solution for Tij in 
various situations ranging from tunneling to ballistic transport. We summa- 
rize the procedure for the linear response calculation. First, we calculate the 
complete Green's function in a tight-binding model. Since the Hamiltonian 
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fJ^i H2 Hi Hi fiN 



Fig. 4. Pictorial representation of the D Amato-Pastawski model for the case of a 
linear chain. 

is sparse this is relatively inexpensive if one uses the decimation techniques 
described in Ref. [18]. With the Green's functions, we evaluate the transmit- 
tances between every pair of sites in the sample (i.e. nodes in the discrete 
equation) and write the transmittance matrix W. Then, we solve for the cur- 
rent conservation equations that involve the inversion of W. 

What are the hmitations of this model? A conceptual one is the momentum 
demolition produced by the locahzed scattering model. Therefore, by decreas- 
ing T^, the dynamics is transformed continuously from quantum ballistic to 
classical diffusive. To describe the transition from quantum ballistic to clas- 
sical ballistic, one should modify the model to have the scattering defined in 
phase space or energy basis. While the first is well suited for scattering matrix 
models [25] , the last is quite straightforward as will be shown in next Section. 
The other aspect is merely computational. Since the resulting matrix W is no 
longer sparse, this inversion is done at the full computational cost. A physi- 
cally appealing alternative to matrix inversion was proposed in DP. The idea 
was to expand the inverse matrix in series in the dephasing collisions, resulting 
in: 

Tr,l = Tk,l + Yl TR,igiTi^L + TK,igiTi,jgjTi^h (36) 

i i j 

+ TR,i9iTi,jgjTi,igiTi;i, + . . . 

i j I 

The formal equivalence with the self-energy expansion in terms of locators 
or local Green's function justifies identifying locator for the classical 

Markovian equation for a density excitation [26] generated by the transition 
probabilities T's. Notice that Eq. (36) can also be rearranged as: 



N ^ 

^R,L = 3r_L + X] '^^,i9iTi,L- (37) 

i=l 

This has the structure of the Dyson equation, graphically represented in Fig. 
5. We notice that according to the optical theorem gi = Ttf,2nhNi, while both 
transmittances entering the vertex are proportional to l/r^, the whole vertex 
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L • = 



R» = 



Fig. 5. Feynman Diagram for the Dyson Equation of the transmittance. It is equiv- 
alent to a particle-hole Green's function in the ladder approximation where a rung 
is represented by a dot. 

is proportional to the dephasing rate. The arrows make explicit that trans- 
mittances are the product of a retarded (electron) and an advanced (hole) 
Green's function. Obviously, one can sum the terms of Eq. (36) to obtain the 
result of Eq.(6) for phenomenological decoherence. 

Many of the results contained in the D'Amato-Pastawski paper for ordered 
and disordered systems were extended in great detail in a series of papers by 
S. Datta and collaborators and are presented in a didactic layout in a book 
[27]. In the next section, we will illustrate how the previous ideas work by 
considering again our reference toy model for resonant tunneling. 



3.3 Effects of decoherence in resonant tunneling 



After an appropriate decimation [18] the "sample" is represented by a single 
state [15,30]. If we choose to absorb the energy shifts into the site energies 
E0 — E0 + A, the Green's function is trivial 



- e-Eo + i{^r+ Rr+ *r)- ^^^^ 

By taking the F's independent on e in the range of interest we get the "broad- 
band" limit. From now on we drop unneeded indices and arguments. From 
this Green's function all the transmission coefficients can be evaluated at the 
Fermi energy. 

(39) 



2r,l 


= 4^r 


K^o,o|^ 


^F, 






1^0,0^ 


^F, and 




= 4^r 


|Go,o|^ 


4>-p 
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We obtain the energy dependent total transmittance: 



The first term in the curly bracket is the coherent contribution while the second 
is the incoherent one. We notice that the effect of the decoherence processes 
is to lower the value of the resonance from its original one in a factor: 



f^r +^ r) 

(41) 

(Lr+ Rr+ -^r) ^ ^ 

In compensation, transmission at the resonance tails becomes increased. 

It is interesting to note the results in the non-linear regime when the voltage 
drop eV is greater than the resonance width (^F + + '^F). If the new 
resonant level lies between eV and yUo, we can easily compute the non- 
linear current using T^x{e,e\/). Notably, one gets that the total current does 
not change as compared with that in absence of decoherent processes, i.e.: 



I- = / rR,L(£, eV)d£ = / Tl^ie, e\/)de (42) 
= 47r^r^j^^,^^r. (43) 

Thus, in this extreme quantum regime, the decoherence processes do not affect 
the overall current. Such relative "stability" against decoherence is fundamen- 
tal in the Integer Quantum Hall Effect [24] and should also be present for 
tunneling through molecular states as both of them have a discrete spectrum. 
Notice that in this last case the dependence of the escape rates on eV is gener- 
ally weak. Hence, the experimental value of the current allows an estimation 
of the escape rate. 

We learn important lessons from the case of resonant tunneling with the in- 
clusion of external degrees of freedom as decoherence: 

1) The integrated intensity of the elastic (coherent) peak is decreased. 

2) An inelastic current comes out to compensate this loss and maintains the 
value of the total transmittance integrated over energy. 

3) Decoherence broadens both contributions to the resonance, relaxing energy 
conservation. 
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The representation of the electron-phonon through complex self energies damps 
the quantum interferences associated with repeated interactions with some vi- 
brational modes that would originate the polaronic states. In what follows 
we will explore some simple Hamiltonian models where decoherence is not 
introduced in such an early stage. Instead of calculating transition rates, the 
approach will be to consider the many-body problem and to compute the 
quantum amplitudes for each state in the Fock space. 



4 The electron-phonon models 

4.1 A state conserving interaction model 



Let us consider a simple model that complements that of DP by providing an 
explicit description of a single extended vibrational mode whose quanta in a 
solid would be optical phonons, 



Hph = hjoh-^h. (44) 

The electrons are described by: 



Tie = E EiCt^i - E [y^,3nh + ^i.^S+cJ . (45) 
i=-oo 

The orbitals between 1 and N = L^/a define our region of interest. Orbitals at 
the edge are connected with the electrodes, where Vij = V, through tunneling 
matrix elements Vo,i = Vi^o = and Vjv.iv+i = Vjv+i.at = ^r- The results 
will be simpler when Vl(r) -C V. Electrons and phonons are assumed to be 
coupled through a local interaction: 



N 



7ie-p/. = E-n^c,(6+ + 6). (46) 



i=l 



To build the electron-phonon Fock space we consider a single electron propa- 
gating in the leads while the number n of vibronic excitations is well defined. 
While this election neglects the phonon mediated electron-electron interac- 
tion, it still has non-trivial elements that are the basis for the development 
of the concept of phonon laser (SASER) [31]. This model will be called State 
Conserving Interaction (SCI). 
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Fig. 6. Scheme for the problem of an electron plus a single plionon mode. The in- 
teracting system in the upper figure is transformed into polaronic modes associated 
to spatial eigenstates which interact weakly through the leads (lower figure). 

The total energy, which is conserved during the transport process, is 



E = en + nkujo with < < ep- 

Here £„ is the kinetic energy of the incoming electron on the left lead and n is 
the number of phonons before the scattering process. A scheme of the complete 
electron-phonon Fock space is shown in Fig. 6. It is clear that an electron that 
impinges from the left when the well has n phonons, can escape through the left 
or right electrodes leaving behind a different number of phonons. Since these 
are physically different situations, the outgoing channels with different number 
of phonons are orthogonal and therefore cannot interfere. This is represented 
in the fact that the sum of transmittances over final channels satisfy unitarity. 

The physical analysis of the excitations can be simplified resorting to a new 

basis to refer the electron-phonon states. We notice that, if we neglect the in- 
teraction with the leads, we can diagonalize the electronic Hamiltonian with- 
out affecting the form of the electron-phonon interaction. First we diagonalize 
the electronic system finding the annihilation operators = 12,f=i Ua,iCi at 
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eigenstates ■^q, with energies Eo. and wave functions Ua,i — {i\ a) 



TV ^ ^ 

no=T. EaK.Ca and He-ph = E -^52^^,(6+ + h) (47) 

a=l 1=0 

The Hamiltonian Hg-ph is just a hnear field for the harmonic oscillator where 
the field amplitude is proportional to the density of the electrons that dis- 
turbs the lattice. Hence, the excitations, called Holstein polarons, are easily 
obtained. The interesting point is that with the proposed Hamiltonian one 
can diagonalize simultaneously every electron subspace. i.e. in this model the 
phonons do not cause transitions between electronic levels. The new excita- 
tions are described by 

OO / 

^k) K,k) ^ Xn,n' {b+) c+|0) (48) 
n'=0 

which is valid at every electron space index a. These operators represent po- 
larons in the energy basis analogous to the Holstein's local polaron model 
[28]. The "polaronic" ground state is related to the unperturbed state by 
|Op,fc) = Xo,n (^^) |0) . The polaronic energies are lower than those of 
electrons plus phonons: 

Ea,n' ^ E^ + fuVo{n' + ^) ~ with n' = (a+a) . (49) 

Since phonons do not produce transitions between electron energy states, this 
model introduces decoherence through a State- Conserving Interaction (SCI). 
The lesson we learn from this is that by writing the interactions in differ- 
ent basis we can choose the quantum numbers that are conserved: e.g. local 
densities in the DP model and energy eigenstates in the SCI model. 

4-2 Coherent and decoherent effects in electronic transport. 

The general problem of the electronic transport has been solved within the 
formalism of Keldysh [15,30] in a FOR decoherent approximation [32]. How- 
ever, here we will solve the full many-body problem. Figure 6 makes evident 
that the phonon emission/absorption can be viewed as a "vertical" hopping in 
a two dimensional network. Once we notice that the Fock-space is equivalent 
to the electron tight-binding model with an expanded dimensionality [31] [33], 
we see that the transmittances can be calculated exactly from the Schrodinger 
equation. While excitations are best discussed in the polaronic basis, for the 
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Fig. 7. Each site is a state in the Fock space: The lower row represents electronic 
states in different sites with no phonons in the well, the sites in black are in the 
leads. Higher rows correspond to higher number of phonons. Horizontal lines are 
hoppings and vertical lines are e-ph couplings. 

electron transport it is preferable the use of the asymptotic states. There, when 
the charge is outside the interacting region, the electron-phonon product states 
constitute the natural basis. Recently, this model has gained additional inter- 
est as it has been used to explain inelastic effects in STM through molecules 
[34] , to study the transport in molecular wires [35] and to investigate Peierls's 
like distortions induced by current in organic compounds [36] . 

A transport calculation is simpler if we prune the Fock space to include only 
states within some range of n allowing a non-perturbative calculation which 
can be considered variational in n. Thus, we are not restricted to a weak e- 
ph coupling. To obtain the transmittances between different channels several 
methods can be adopted. One possibility is to solve for the wave function 
iteratively [33]. An alternative is to obtain Green's functions to get the trans- 
mittances. In this case, the horizontal "dangling chains" in the Fock space 
can be eliminated through a decimation procedure [12,18] introducing com- 
plex self-energies in the corresponding "sites". 

For simplicity, let us consider the case of a single state of energy Eq in the re- 
gion of interest which we will call the "resonant" state. It could be interpreted 
as a HOMO or LUMO state depending on the situation. It interacts with a 
dispersionless phonon mode, and it is coupled to a source and drain of charge. 
The problem for an electron that tunnels through the system can be mapped 
to the one-body problem shown in Fig. 7. 

Let us consider an asymptotic incoming scattering state consisting of a wave 
packet built with electron-phonon states in the left branch corresponding to 
Uq phonons, i.e. an electron coming from the left while there are tiq phonons in 
the well. When it arrives at the resonant site where it couples to the phonon 
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field, it can either keep the available energy E as kinetic energy or change it 
by emitting or absorbing n phonons. Each of these processes contributes to 
the total transmittance which is given by: 

n 

In Fig. 8-a we show the total transmittance (thick line) for where is 

zero and the e-ph coupling is weak. There, the appearance of satellite peaks 
at energies Eo + nhwo can be appreciated. To discriminate the processes con- 
tributing to the current, we include with a dashed line, the elastic contribution 
to the transmittance, i.e. that due to electrons which are escaping to the right 
without leaving vibrational excitations behind. This makes almost all of the 
main peak and just decreasing portions of the satellite peaks. The elastic 
contribution at the satellite peaks corresponds to a sum of virtual processes 
consisting of the emission of phonons followed by their immediate reabsorp- 
tion. The inelastic component associated with the emission of one phonon is 
shown with a dotted line. 

While the total transmittance (thick line) shows a smooth behavior, the elas- 
tic component (dashed line) exhibits a strong dip in the region between the 
first two resonances. Therefore, almost all of the transmitted electrons within 
this energy range will be scattered to the inelastic channels. This sharp drop 
in the elastic transmittance is produced by a destructive interference between 
the different possible "paths" in the Fock-space connecting the initial and the 
final channels. These paths can be classified essentially as a direct "path" be- 
tween the incoming and outgoing channels and the same path dressed with 
virtual emission and absorption processes. The first inelastic component (dot- 
ted line) also shows a similar behavior in the valley between the second and 
third resonances. The main factors that control the magnitude of these antires- 
onances are the escapes to the leads and the e-ph coupling. This concept was 
introduced in Refs. [37,18] to extend the Fano-resonances [38] observed in spec- 
troscopy to the problem of conductance. Perfect antiresonances correspond to 
situations where both the real and the imaginary part of the transmission 
amplitude are zero. In the present case, such a perfect interference condition 
does not occur, and there is a non-zero minimum transmission at the dips. 
Through a similar analysis one could tailor the geometrical parameters of the 
system (as done in Ref. [39]) to optimize the phonon emission. 

An alternative way to appreciate this interference effect is to plot the path 
followed by the transmission amplitude in the complex plane when the ki- 
netic energy of the incoming electron is changed [40,41]. A plot for the elastic 
transmission amplitude is shown in Fig. 8-c). For e = 0, there is no trans- 
mission. If the energy is increased, one starts to follow the path in the figure 
anti-clockwise. After reaching the point corresponding to the first resonance, 
the transmission starts to decrease and the curve develops a turning near the 
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Fig. 8. a) Transmittance as a function of the incident electronic kinetic energy. The 
total transmittance is shown with a thick solid line, the elastic component with a 
dashed line and Ti^o with a dotted line. The non-interacting transmittance is also 
shown with a thin line. In b) the phase shift of the elastic transmittance as a function 
of the energy and for the non-interacting case (dotted line) are shown. Figure c) 
shows the path of the elastic transmission amplitude in the complex plane when 
the electronic energy is changed. The full circles correspond to the first and second 
peaks shown in a). These results are obtained with Eq = —1.5, Vj, = Vr = 0.1, 
Vg = 0.1 and Huq = 0.2. 

origin. The first antiresonance takes place at the point of minimum distance 
from the origin. 



In order to rationalize the main processes involved in the first two peaks, let 
us represent them schematically in Fig. 9. Panel (a) shows the standard elastic 
process in which no phonon is emitted. Panel (b) is a notable effect that occurs 
when the phonon emission is virtual. Notice that this virtual process can have 
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Fig. 9. Schematic representation of the processes that lead to the elastic and inelastic 
components of the peaks in the total transmittance. Figures a) and b) correspond 
to the first and second peaks of the elastic transmittance respectively. Figures c) 
and d) correspond to the inelastic part. The well's ground state is represented by 
a solid line and the first excited polaron state by a dotted line. The final polaron 
states at the right of the well are also shown. These states have an energy equal to 
the incident electron kinetic energy. The final polaron state is depicted with a solid 
line for the elastic case and with a dotted line for the inelastic situation. The level 
corresponding to the electron final energy is also shown as a solid line for this case. 

strong consequences. In Fig. 8 it produces an increase of the transmittance 
in almost two orders of magnitude. Panel (d) shows a real inelastic process 
which is expected to give a peak if the initial energy satisfies E = Eo + huo- 
Notably, one must expect also a contribution when e = Eg which corresponds 
to the virtual tunneling into the resonant state before emitting a phonon as 
shown in panel (c). 



Another interesting quantity that we can explore with the present formalism 
is the phase of the transmitted electron through different channels, 



1 

2i G m 



(51) 



whose energy derivative gives information on the dynamics of the process. 
Figure 8-b shows the phase shift of the elastic transmission probability. For 
reference, the phase shift in the absence of e-ph interaction is also shown 
with a dotted line. It increases by vr over the width of the non-interacting 
transmission resonance. The same occurs in the interacting case. However, it 
can be seen that each satellite peak has associated a phase fluctuation. Instead 
of an increase in the phase by tt which would occur for real resonant peaks, 
across each "satellite peak" associated with the virtual processes, there is a 
phase dip that results in consecutive resonances having phases close to -7r/2. 
These phase dips are a manifestation of the anti-resonances shown in Fig. (8- 
a) for the elastic transmittance. For perfect zero transmission points, one has 
an abrupt phase fall of vr instead of the smooth phase slip shown in the Fig. 
(8-b) [40]. 
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Let us consider a more general case in which there are initially tiq phonons in 
the scattering region. The vertical hopping matrix element connecting states 
with hq and ng + 1 phonons is \/no + IVg. Then, under these conditions, the 
adimensional parameter. 



characterizes the strength of the e-ph interaction. It presents two regimes 
according to the importance of this interaction. 

In the limit ^ <^ 1 the "vertical" processes are in a perturbative regime. The 
"vertical hopping" y/rio + IVg will not be able to delocalize the initial state 
along this "vertical direction" and therefore the elastic contribution is the 
most important. In Fig. (10-a), we show with a thick line the transmission 
probability for a case where g — 0.25 in presence of 10 phonons. Notice that 
"satellite" resonances corresponding to phonon emission and absorption are 
separated by hujo from the main resonance. For comparison, the curve in ab- 
sence of elcctron-phonon interaction is also included with a dotted line. We see 
that, as with the DP model, the main resonance peak decreases and presents 
a general broadening where we can recognize details of the excitation struc- 
ture of the phonon field. The phase shift as a function of the electronic energy 
is shown in b). The solid line is the phase shift for the elastic component of 
the total transmittance. There, we can appreciate the phase fluctuations that 
appear even in the elastic transmission. 

On the other hand, if ^ ^ 1, the e-ph interaction is in the non-perturbative 
regime. This can be achieved either by a strong Vg, or by a high n. The 
second situation would correspond to the case of high temperatures or a far 
from equilibrium phonon population. In that case, the total transmittance can 
be obtained by summing the transmittances for the different possible initial 
conditions weighted by the appropriate thermal factor [42]. In any of these 
situations, the resulting strong e-ph coupling leads to a breakdown of pertur- 
bation theory. A similar case is found for the quantum dots studied in Ref. 
[43]. There, in contrast to the situation in bulk material, the studied quan- 
tum dots show a strong e-ph coupling. In molecular wires, where elementary 
excitations (e or ph) are highly confined, one can expect a similar scenario. 

In Fig. (11-a), we show the transmittance as a function of the incident elec- 
tron kinetic energy for the case where Vg is strong 4 ) and there are no 
phonons in the well before the scattering process. In contrast to the perturba- 
tive regime, where the energy shift of the peaks {—Vg/huJo = —ghujo) is small, 
here the energy shift is important and the inelastic contributions to the total 
transmittance dominate. In the limit of weak coupling between the scattering 
region and the leads, the transmission probability through the channel with n 




(52) 
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Fig. 10. a) Total transmittance as a function of the incident electronic kinetic energy. 
The solid line corresponds to the case in which there are 10 phonons in the well 
before the scattering process. The dotted line corresponds to the case in which there 
is no interaction with the vibrational degrees of freedom, b) Phase shift in units of 
TT as a function of the incident electronic kinetic energy. The solid line is the phase 
shift for the elastic transmission. The dotted line corresponds to the non-interacting 
case. The parameters of the Hamiltonian in units of the hopping V are: Eq = 0, 
Vl = Vr = 0.05, Vg = 0.004 and hwo = 0.02. 

phonons is: ex.p{—g)g"'/n\. 

Figure (H-b) shows the total transmittance as a function of the incident elec- 
tron kinetic energy for an extreme case where there are 300 phonons in the 
well before the scattering process. The electron can absorb or emit as many 
phonons, N^jf = A^ynVg/huj^j, as allowed by the interaction strength. This 
means that the phonon spectrum, weighted on the local electron-phonon state 
n, will be quite independent on the details of the spectral densities at the far 
end of the effective "vertical" chain i.e. the states n±Neff. Hence, the spectral 
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Fig. 11. a) and b) show the total transmittance as a function of the incident kinetic 
electronic energy. The area under the elastic transmittance is shaded in gray. As 
before, the dotted line corresponds to the non- interacting case. In a), the solid 
line corresponds to the case in which there are no phonons in the well before the 
scattering process and Vg = 0.100, Huq = 0.05. In b) , the solid line corresponds to 
the case in which there are 300 phonons in the well before the scattering process 
and Vg = 0.0015 and Huiq = 5.010"^. The parameters in units of the hopping V are: 
Eo = 0,Vl = Vr = 0.05. 

density will show the typical form of a one a dimensional band even when we 
take Ngff = 4:(yg/hwo)y/n <C N. This feature manifests in the transmittance 
presented in Fig. 11-b). The unperturbed resonant peak is a dotted line which 
contrasts with the total transmission in presence of phonons shown with a 
continuous line. Its trace follows the structure of the phonon excitation with 
the typical square root divergences at the band edges smoothed out by the 
inhomogeneity in the hopping elements and the uncertainty introduced by the 
escape to the leads. 
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5 The solution of Time Dependence 

Finally, we would like to present the basic features of time dependent trans- 
port. This is relevant since a coherent "time of flight" should be shorter than 
any conformational correlation time [44] . The basic idea to get the physics of 
time dependent phenomena is to obtain the evolution of an arbitrary initial 
boundary condition of the form ['4'{Xj)i'*{Xk)] source ■ Since here Xj = {rj,tj) 
with general tj, this essentially generalizes a density matrix which introduces 
temporal correlations that define the energy of this injected particle. The 
"density" at a later is obtained from the exact solution of the Schrodinger 
equation: 



(53) 

In order to establish a correspondence with the Danielewicz solution to the 
Schrodinger equation in the Keldysh formalism, we used the continuous vari- 
able Green's function which is related to the exact discrete one in the open 
system by 

G^Kr^.e) ^^G^^i{e)Mri)vKrk). (54) 

k,l 

Now, the key is to recognize that in any Green's function, a macroscopically 
observable time is t — ^[tj + tk] (time center). Its Fourier transform is an 
observable frequency u. Meanwhile, energies are associated with internal time 
differences tj —tk (time chords). Within this scheme, the time correlated initial 
condition (source), determining the occupation of a local orbital is expressed 
in terms of the time independent local density of states Ni{e) : 

[ui(tK(t)Lurce = J^l^^t [u,{t + 5t/2)u:(t-5t/2)l^^^^^eMis5t] 

= j deNi{e)i{e,t), (55) 

and the occupation factor f(£, t). Introducing this notation in Eq.(53) one can 
identify the dynamical transmittances Ti j[£,Lo). We refer to [15] for a detailed 
manipulation of the time integrals. The basic result is 



r,,-(£,u;) = 2r,(£)G£.(£ + \nw)2T,{e)G%{e - \hw). (56) 
The time dependent transmittance is then: 
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(57) 



where one recovers our old steady state transmittance as: 




Equation (53) then becomes 




which is the Generahzed Landauer-Biittiker Equation (GLBE) [26,15]. Ac- 
cording to Eq.(58), the first term accounts for the particles that are leaving 
the reservoir at site j at time t to suffer a dephasing collision at some future 
time at orbital i. The second accounts for particles that having had a previous 
dephasing collision at time ti at the i orbital reach site j at time t. 

In the GLBE formulation the essential features of time dependence in trans- 
port is contained in the transmittances. Since the spectrum is continuous, we 
can keep the lowest order in the frequency expansion e.g. 



According to Eq. (4.5) in Ref. [15] the propagation time rp is identified with 
the first significant coefficient in this expansion. Typically, it results the first 
order. 



which can be identified with the Wigner time delay. This propagation time was 
evaluated in various simple systems in Ref. [15] recovering the ballistic and 
diffusive times for clean and impure metals respectively. In a double barrier 
system, in the resonant tunneling regime, the propagation time is determined 
by the life-time inside the well. In fact, using the functions of Subsection 2.3, 
one gets for the propagation through the resonant state: 




(60) 



-1 



(61) 



2 de [Gf.ie) 



(62) 



Tp = 



2(^r r) ■ 
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Prom these considerations, we see that Tp represents a hmit to the response in 
frequency (admittance) of the device. Typically, one gets = Go/(l — iouTp). 
This is in fair agreement with the experimental results [45]. 

As an striking example, we mention the "simple" case of tunneling through a 
barrier [46] of length L and height U exceeding the kinetic energy e of the 
particle. For barriers long enough the expantion is dominated by the second 
order term in Eq. (60) and one gets: 



which, within our non-relativistic description, can be extremely short provided 
that the barrier is high enough. This is the time that one has to compare with 
vibronic and configurational frequencies [47]. 

The general propagation times can also be calculated in more complex situa- 
tions such as disordered systems [48] and those affected by incoherent inter- 
actions [15]. 



6 Perspectives 

We have presented the general features of quantum transport in mesoscopic 
systems. There, interactions with environmental degrees of freedom introduce 
complex phenomena whose overall effect is to decrease the clean interfer- 
ences expected for an isolated sample. Such degradation can be accounted 
for as decoherence. We introduced various simple models which present de- 
coherence. We started presenting the simplest phenomenological models of 
Biittiker, by discussing its connection with the Hamiltonian description of the 
DP model. Various simple models for resonant tunneling devices including 
electron-phonon interactions were introduced. We feel that they contain the 
essential coherent and decoherent effects induced by the vibronic degrees of 
freedom. 

Reference [12] showed that in disordered systems, the effective transmittance 
away from the resonances keeps its form as a superposition of tails from dif- 
ferent resonances. This applies to both the DP and SCI models. This justifies 
the simplification of using a single resonant state. 

Our results for this model make quite clear how complex many body inter- 
actions result in the loss of the simple interferences of one body description. 
Essentially, each external degree of freedom coupled with the electronic states 
leads to two situations producing decoherence: 




(63) 
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1) By real emission or absorption of phonons, they open an additional scat- 
tering channels contributing incoherently to the transport. This situation is 
represented in the DP decoherence model, as well as in the various polaronic 
models presented here. Those are real processes which can be detected by 
measuring a change in the bath state (phonon models) or energy dissipated 
(DP model) within the sample. 

2) However, even when the electron-phonon processes might be virtual , they 
would add new alternatives to the quantum phase of the outgoing state. This 
strongly limits an attempt to control the electron phase and hence is mani- 
fested as decoherence. 

The study of the coherent component in the presence of decoherent processes 
gives a first, though imperfect, hint to the conditions of the transport pro- 
cesses. In particular, it can show how decoherence can affect the interference 
between different propagation pathways. Once they are summed up, the con- 
ductance, which is a square modulus, would manifest the effect of a diminished 
interference. This approach was adopted in a recent work [49] that addressed 
the delocalizing effect of decoherence on transport in the variable range hop- 
ping regime. 

It is also worthwhile to mention that the randomization of the quantum phase 
introduced by the virtual processes can be as much effective [50] as those 
involving a real energy exchange in producing decoherence. In fact, it has 
been found recently that this practical uncertainty is the mechanism through 
which quantum chaos contributes to dissipation and irreversibility [51]. Such 
quantum chaotic systems have intrinsic decoherence time scales which contrast 
with the extreme quantum regime of a resonant tunneling where decoherence, 
if any, is controlled by the environment. 

We want to close this section mentioning the connections of the results in Sec- 
tion 4 with other ongoing research on various hot issues of electronic transport. 
While the problem of tunneling times [52,46] is a controversial one, it has im- 
portant practical aspects [45]. Indeed, molecular electronics opens the whole 
issue of quantum dynamical processes to a fresh consideration. One aspect is 
the effect of decoherence on frequency response. A related issue under study is 
the interconnection between decoherence times, irreversibility and dynamical 
chaos [50]. Having shown the subtle relation between spectral properties and 
time dependences, one foresees that our models can provide new insight to this 
topic. The analysis of the consequences of the spectral complexity of many- 
body systems is a fully unexplored field ahead. Once again, technology pushes 
us to the frontiers of the conceptual understanding of Quantum Mechanics 
and Statistical Physics. 



29 



7 Acknowledgments 



We acknowledge P. R. Levstein and F. Toscano for critical discussions. We 

acknowledge financial support of bi-national programs of Antorchas-Vitae and 
CONICET-CONICIT as well as local support from SeCyT-UNC. HMP and 
LEFFT are affiliated with CONICET. 



References 

[I] C. Joachim, J. K. Gimzcwski and A. Aviram. Nature, 408 (2000) 541; M.A. 
Reed and J. M. Tour, Sci. Am., 282 (6) (2000) 68. 

[2] C. Dekker and M. Ratner, Physics World 14 (2001) 29. 

[3] H. Prauenfelder, R G. Wolynes and R. H. Austin, Rev. Mod. Phys. 71 (1999) 
S419. 

[4] P. W. Anderson, Rev. Mod. Phys. 50 (1978) 191 

[5] L. Esaki, Rev. Mod. Phys. 46 (1974) 237. 

[6] M. Kastner, Physics Today 46 (1993) 25. 

[7] S. Washburn and R. A. Webb, Rep. Prog. Phys. 55 (1993) 1311. 

[8] V. J. Goldman and D. C. Tsui and J. E. Cunningham, Phys. Rev. B, 36 (1987) 
7635. 

[9] B. C. Stipe, M. A. Rezaei, and W. Ho, Phys. Rev. Lett. 81 (1998) 1263 ; H. 
Park et al., Nature 407 (2000) 57. 

[10] J. H. Schon, A. Dodabalapur, Z. Bao, Ch. Kloc, O. Schenker, B. Batlogg, Nature 
410, (2001) 189. 

[II] R. Landauer, Phys. Lett. 85A (1981) 91; R. Landauer, IBM J. Res. Develop. 
1 (1957) 223 ; R. Landauer, Philos. Mag. 21 (1970) 863. 

[12] J. L. D'Amato and H. M. Pastawski, Phys. Rev. B 41 (1990) 7411. 

[13] M. Biittiker, Phys. Rev. Lett. 57 (1986) 1761. 

[14] V. L. Nguyen, B. Z. Spivak, and B. L Shklovskh, Pisma Zh. Eksp.Teor. Fiz. 
41(1985) 35 [JETP Lett. 41 (1985) 42]; E. Medina and M. Kardar, Phys. Rev. 
B 46 (1992) 9984. 

[15] H. M. Pastawski, Phys. Rev. B, 46 (1992) 4053. 

[16] M. Biittiker. IBM J. Res. Develop. 32 (1988) 63. 



30 



[17] P. O. Lowdin, J. Chem. Phys. 19, 1396 (1951); these projection techniques were 
discussed in the context of a real space Renormalization Group Decimation in: 
C.E.T. Gongalves da Silva and B. Koiller, Sohd St. Comm. 40 (1981) 215 and 
J. V. Jose, Phys. Rev. Lett. 49 (1982) 334. 



[is: 

[19 

[2o; 

[21 
[22; 

[23; 

[24 
[25; 
[26 
[27 

[28 
[29 
[30 

[31 

[32; 
[33 

[34 
[35; 



P. Levstein, H. M. Pastawski and J. L. D'Amato, J. Phys.: Condens. Matter 2 

(1990) 1781. 

H. M. Pastawski and E. Medina, Rev Mex. Fis. 47S1 (2001) 1; also available 
at |cond-mat/0103219| . 

L. V. Keldysh, Zh. Eksp. Teor. Fiz. 91, (1986) 1815 [Sov. Phys. JETP 64 
1075 (1986)]; L. P. KadanofF and G. Baym, Quantum Statistical Mechanics 
(Benjamin, New York, 1962). 

P. Danielewicz, Ann. Phys. (N. Y.) 152 (1984) 239. 

D. S. Fisher and P. A. Lee, Phys. Rev. B 23, 6951 (1981); F. Sols, Ann. Phys. 
(N.Y.) 214 (1992) 386. 

J. N. Onuchic, P.C. P. de Andrade and D. N. Beratan, J. Chem. Phys. 95 (1991) 
1131; V. Mujica, M. Kemp and M.A. Ratner, J. Chem. Phys. 101 (1994) 6856; 
M.D. CoutinhoNeto and A.A. daGama, Chem. Phys. 203, 43 (1996). 

F. Gagel and K. Maschke, Phys. Rev. B 54 (1996) 13889. 

I. Knittel, F. Gagel and M. Schreiber. Phys. Rev. B, 60 (1999) 916. 
H. M. Pastawski Phys. Rev. B 44 (1991) 6329. 

S. Datta, J. Phys. Cond. Matter 2 (1990) 8023; Electronic Transport in 
Mesoscopic Systems, Cambridge U. P. (Cambridge, 1995) and references 
therein. 

N. S. Wingreen, K. W. Jacobsen, and J. W. Wilkins, Phys. Rev. Lett. 61 (1988) 
1396; Phys. Rev. B, 40 (1989) 11834. 

J. A. St0vneng, E. H. Hauge, P. Lipavsky and V. Spicka, Phys. Rev. B 44 

(1991) 13595. 

H. Haug and A. -P. Jauho. Quantum Kinetis in Transport and Optics 
of Semiconductors. Springer- Verlag (Heidelberg, 1998); A. P. Jauho, N.S. 
Wingreen and Y. Meir, Phys. Rev. B 50 (1994) 5528. 

E. V. Anda, S. S. Makler, H. M. Pastawski, and R. G. Barrera, Braz. J. Phys. 
24 (1994) 330. 

R. G. Lake, G. Khmeck and S. Datta, Phys. Rev. B 47 (1993) 6427. 

J. Bonca and S. A.Trugman, Phys. Rev. Lett, 75, 2566 (1995); Phys. Rev. Lett. 
79 (1997) 4874. 

N. Mingo and K. Makoshi, Phys. Rev. Lett. 84 (2000) 3694. 
E. G. Emberly and G. Kirczenow, Phys. Rev. B 61 (2000) 5740. 



31 



[36; 
[37; 
[38; 

[39 

[4o; 

[41 

[42; 
[43; 



[44; 
[45; 

[46; 

[47; 

[48; 



[49; 
[50; 

[51 



H. Ness, S.A. Shevlin and A.J. Fisher, Phys. Rev. B 63 (2001) 125422. 

J. L. D'Amato, H.M. Pastawski and J. F. Weisz, Phys. Rev. B 39 (1989) 3554. 

U. Fano, Phys. Rev. 124 (1961) 1866. 

L. E. F. Foa Torres, H. M. Pastawski, and S. S. Makler, Phys. Rev. B. 64 (2001) 
193304. 

H. W. Lee, Phys. Rev. Lett. 82 (1999) 2358. 

T. Taniguchi and M. Biittiker, Phys. Rev. B 60 (1999) 13814. 

K. Haule and J. Bonca, Phys. Rev. B 59 (1999) 13087. 

S. Hameau, Y. Guldner, O. Verzelen, R. Ferreira, G. Bastard, J. Zeman, A. 
Lemaitre and J.M. Gerard, Phys. Rev. Lett. 83 (1999) 4152; O. Verzelen, S. 
Hameau, Y. Guldner, J.M. Gerard, R. Ferreira, and G. Bastard, Jpn. J. Appl. 
Phys. 40 (2001) 1. 

M. A. Ratner, PNAS, 98 (2001) 387. 

E. R. Brown, T. C. L. G. Sollner, C. D. Parker, W. D. Goodhue, and C. L. 
Chen, Appl. Phys. Lett. 55 (1989) 1777. 

M. Biittiker and R. Landauer, Phys. Rev. Lett. 49 1739 (1982); Physica Scripta, 
32 (1985) 429; R. Landauer and Th. Martin, Rev. Mod. Phys. 66 (1994) 217. 

A. Nitzan, J. Jortner, J. Wilkie, A. L. Burin and M. Ratner, J. Phys. Chem. 

B, 104 (2000) 5661. 

Prigodin, V. N. , B. L. Altshuler, K. B. Efetov and S. lida, Phys. Rev. Lett. 72 
(1994) 546; Pastawski, H. M., P. R. Levstein and G. Usaj, Phys. Rev. Lett 75, 
4310 (1995); Pastawski, H. M., G. Usaj and P. R. Levstein, Chem. Phys. Lett. 
261 (1996) 329. 

E. Medina and H. M. Pastawski, Phys. Rev. B 61 (2000) 5850. 

R. A. Jalabert, H. M. Pastawski, Phys. Rev. Lett. 86 (2001) 2490. 

F. M. Cucchietti, H. M. Pastawski and R. Jalabert, Physica A 283 (2000) 285; 
F. M. Cucchietti, H. M. Pastawski and D. A. Wisniacki, to appear in Phys. 
Rev. E (also available at jond-mat /0102135 ) . 



[52] C. J. Bolton-Heaton, C. J. Lambert, V. I. Falco, V. Prigodin and A.J. Epstein, 
Phys. Rev. B 60 (1999) 10569. 



32 



